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Chaos analyses in both phases of QED and QCD 
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We analyze the leading Lyapunov exponents of U(l) and SU(2) gauge field configurations on the lattice which 
are initialized by quantum Monte Carlo simulations. We find that configurations in the strong coupling phase are 
substantially more chaotic than in deconfinement. 



1. Motivation 

The study of chaotic dynamics of classical field 
configurations in field theory finds its motivation 
in phenomenological applications as well as for 
the understanding of basic principles. The role 
of chaotic field dynamics for the confinement of 
quarks is a longstanding question. Here, we an- 
alyze the leading Lyapunov exponents of com- 
pact U(l) and of SU(2)-Yang-Mills field configu- 
rations on the lattice. The real-time evolution of 
the classical field equations was initialized from 
Euclidean equilibrium configurations created by 
quantum Monte Carlo simulations. This way we 
expect to see a coincidence between the strong 
coupling phase and the strength of chaotic be- 
havior in lattice simulations. 

After reviewing essential definitions of the 
physical quantities describing chaos and their 
computation in lattice gauge theory we out- 
line our method for the extraction of starting con- 
figurations of a three dimensional Hamiltonian 
dynamics from four dimensional Euclidean field 
configurations. Our results are then presented 
by showing an example of the exponential diver- 
gence of small initial distances between nearby 
field configurations. It is followed by a detailed 
study of the maximal Lyapunov exponent and av- 
erage plaquette energy as a function of the cou- 
pling strength. 

2. Classical chaotic dynamics 

Chaotic dynamics in general is characterized 
by the spectrum of Lyapunov exponents. These 



exponents, if they are positive, reflect an expo- 
nential divergence of initially adjacent configu- 
rations. In case of symmetries inherent in the 
Hamiltonian of the system there are correspond- 
ing zero values of these exponents. Finally neg- 
ative exponents belong to irrelevant directions 
in the phase space: perturbation components in 
these directions die out exponentially. Pure gauge 
fields on the lattice show a characteristic Lya- 
punov spectrum consisting of one third of each 
kind of exponents This fact reflects the elim- 
ination of longitudinal degrees of freedom of the 
gauge bosons. Assuming this general structure of 
the Lyapunov spectrum we investigate presently 
its magnitude only, namely the maximal value of 
the Lyapunov exponent, Lmax- 

The general definition of the Lyapunov expo- 
nent is based on a distance measure d{t) in phase 
space. 
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In case of conservative dynamics the sum of all 
Lyapunov exponents is zero according to Liou- 
ville's theorem, 



(2) 



We utilize the gauge invariant distance measure 
consisting of the local differences of energy den- 
sities between two field configurations on the lat- 
tice: 
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Here the symbol stands for the sum over all 
Np plaquettes, so this distance is bound in the 
interval (0, 2N) for the group SU(N). Up and U'p 
are the familiar plaquette variables, constructed 
from the basic link variables U^^i, 

U,,, = exp {aAl^^T'^) , (4) 

located on lattice links pointing from the position 
X = {xi, X2,X3) to X + aci. The generators of the 
group are T"^ = —igT'^/2 with r"^ being the Pauli 
matrices in case of SU(2) and A'^ ^ is the vector 
potential. The elementary plaquette variable is 
constructed for a plaquette with a corner at x 
and lying in the ij-plane as 

Ux,ij = Ux.iUx+ijUl^j^^ulj. (5) 

It is related to the magnetic field strength j,: 



Ux 



exp [e^jkaBl^kT" 



(6) 



The electric field strength E'^ ^ is related to the 
canonically conjugate momentum Px^i = Ux.i via 
2a 
? 



(7) 



3. Initial states prepared by quantum 
Monte Carlo 

The Hamiltonian of the lattice gauge field sys- 
tem can be casted into the form 
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(8) 



Here the scalar product between group elements 
stands for {A,B) = \iv{AB^). The staple vari- 
able y is a sum of triple products of elementary 
link variables closing a plaquette with the cho- 
sen link U . This way the Hamiltonian is formally 
written as a sum over link contributions and V 
plays the role of the classical force acting on the 
link variable U . The naive equations of motion 
following from this Hamiltonian, however, have 
to be completed in order to fulfill the constraints 

([/,[/) = 1, 

{P,U) = 0. (9) 
The following finite time step recursion formula: 
U' = U + dt{P' -eU), 

P' = P + dt{V - fiU + eP'), (10) 



with the Lagrange multipliers 

e = {U,P'), 
fJ- = 



{U,V) + {P',P'), 



(11) 



conserves the Noether charge belonging to the 
Gauss law, 



r = J2PU^-Y.u^p 



(12) 



Here the sums indicated by + run over links start- 
ing from, and those by — ending at a given site 
X, where the Noether charge F is defined. The 
above algorithm is written in an implicit form, 
but it can be casted into explicit steps, so no it- 
eration is necessary Q . 

Initial conditions chosen randomly with a given 
average magnetic energy per plaquette have been 
investigated in past years. In the SU(2) case, a 
linear scaling of the maximal Lyapunov exponent 
with the total energy of the system has been es- 
tablished for different lattice sizes and coupling 
strengths [Q . In the present study we prepare the 
initial field configurations from a standard four 
dimensional Euclidean Monte Carlo program on 
a 12^^ X 4 lattice varying the inverse gauge cou- 
pling /3 |]. 

We relate such four dimensional Euclidean lat- 
tice field configurations to Minkowskian momenta 
and fields for the three dimensional Hamiltonian 
simulation by the following approach: 

First we fix a time slice of the four dimensional 
lattice. We denote the link variables in the three 
dimensional sub-lattice by U' = Ui{x, t). Then we 
build triple products on attached handles in the 
positive time direction, 

U" ^UA{x,t)U^{x,t + a)ul{x + a,t). We obtain 
the canonical variables of the Hamiltonian system 
by using 

P = iU"-U')/dt, 

U oc {U" + U'). (13) 

Finally U is normalized to {U, U) = 1. 

This definition constructs the momenta accord- 
ing to a simple definition of the timelike covariant 
derivative. The multiplication with the link vari- 
ables in time direction can also be viewed as a 
gauge transformation to U4{x,t) — 1, i.e. Aq = 
Hamiltonian gauge. 



3 



4. Chaos and confinement 

We start the presentation of our results with 
a characteristic example of the time evolution of 
the distance between initially adjacent configura- 
tions. An initial state prepared by a standard four 
dimensional Monte Carlo simulation is evolved 
according to the classical Hamiltonian dynamics 
in real time. Afterwards this initial state is ro- 
tated locally by group elements which are chosen 
randomly near to the unity. The time evolution of 
this slightly rotated configuration is then pursued 
and finally the distance between these two evolu- 
tions is calculated at the corresponding times. A 
typical exponential rise of this distance followed 
by a saturation can be inspected in Fig. |l| from 
an example of U(l) gauge theory in the confine- 
ment phase and in the Coulomb phase. While the 
saturation is an artifact of the compact distance 
measure of the lattice, the exponential rise (the 
linear rise of the logarithm) can be used for the 
determination of the leading Lyapunov exponent. 
The naive determination and more sophisticated 
rescaling methods lead to the same result. 

The main result of the present study is the de- 
pendence of the leading Lyapunov exponent Lmax 
on the inverse coupling strength (3, displayed in 
Fig. H for a statistics of 100 independent U(l) 
configurations. As expected the strong coupling 
phase, where confinement of static sources has 
been established many years ago by proving the 
area law behavior for large Wilson loops, is more 
chaotic. The transition reflects the critical cou- 
pling to the Coulomb phase. Furthermore the 
maximal Lyapunov exponent scatters more pro- 
nounced than the average energy per plaquette. 
Fig. ^ shows the somewhat smoother transition 
of the energy per plaquette as a function of the 
inverse coupling strength. Fig. ^ depicts the cor- 
relation of the Lyapunov exponents and the pla- 
quette energies for 100 U(l) configurations. The 
blank area is indicative of the transition point be- 
ing presumable of first order. 

We now turn to a comparison of expectation 
values of U(l) and SU(2) theory. Fig. || exhibits 
the averaged leading Lyapunov exponent between 
the strong and the weak coupling regime. The 
smoother fall-off of the SU(2) Lyapunov expo- 
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Figure 1. Exponentially diverging distance of 
initially adjacent U(l) field configurations on a 
12^ lattice prepared a,t (3 — 0.9 in the confinement 
phase (top) and at /3 = 1.1 in the Coulomb phase 
(bottom). 

nent reflects the second order of the finite tem- 
perature transition to a Debye screened phase of 
free quarks. Fig. ^ compares the averaged pla- 
quette energies of both gauge theories signaling 
the different order of their phase transitions. 

Fig. shows the energy dependence of the Lya- 
punov exponents for both theories. One observes 
an approximately linear relation for the SU(2) 
case while a quadratic relation is suggested for the 
U(l) theory in the weak coupling regime. From 
scaling arguments one expects a functional rela- 
tionship between the Lyapunov exponent and the 
energy 

(14) 



L{a) oc a'=~i£;'=(a), 



with the exponent k being crucial for the contin- 
uum limit of the classical field theory. A value of 
fc < 1 leads to a divergent Lyapunov exponent, 
while A; > 1 yields a vanishing L in the contin- 
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Figure 2. Transition of the leading Lyapunov 
exponents from 100 U(l) configurations as a func- 
tion of the inverse coupling strength (3. 
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Figure 3. Transition of the plaquette energy as 
in Fig. 2. 
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Figure 4. Scatter plots of Lyapunov exponents 
and plaquette energies for 100 U(l) configura- 
tions. 



uum. The case fc = 1 is special leading to a finite 
non-zero Lyapunov exponent. Our analysis of the 
scaling relation (14) gives evidence, that the clas- 
sical compact U(l) lattice gauge theory has A: 2 
and with L(a) — > a regular continuum theory. 
The non-Abelian SU(2) lattice gauge theory sig- 
nals fc « 1 and stays chaotic approaching the con- 
tinuum. 

Summarizing we investigated the classical 
chaotic dynamics of U(l) and SU(2) lattice gauge 
field configurations prepared by quantum Monte 
Carlo simulation. The maximal Lyapunov expo- 
nent shows a pronounced transition as a function 
of the coupling strength. Both for QED and QCD 
we find that configurations in the strong coupling 
phase are substantially more chaotic than in the 
weak coupling regime. Our results demonstrate 
that chaos is present when particles are confined, 
but it persists partly also into the Coulomb and 
quark-gluon-plasma phase. So far the situation 
for the gauge fields. An independent analysis of 
the fermion fields yields compatible results with 
respect to quantum chaos 1^. 
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Figure 5. Comparison of the average maximal Lyapunov exponent in U(l) gauge theory with (3 = 
(left) and in SU(2) gauge theory with (3 = A/g"^ (right) when crossing from the strong to the weak coupling 
phase. 
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Figure 6. Comparison of the average plaquette energy as in Fig. 5. 
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Figure 7. Comparison of average maximal Lyapunov exponents as a function of the scaled average 
energy per plaquette ag'^E. The U(l) theory (left) shows an approximately quadratic behavior in the 
weak coupling regime whereas the SU(2) theory (right) is approximately linear. 



